Nuclear spin cooling using Overhauser field selective coherent population trapping 
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Hyperfine interactions with a nuclear spin environment fundamentally limit the coherence proper- 
ties of confined electron spins in the solid-state. Here, we show that a quantum interference effect in 
optical absorption from two electronic spin states of a solid-state emitter can be used to prepare the 
surrounding environment of nuclear spins in well-defined states, thereby suppressing electronic spin 
dephasing. The evolution of the coupled electron-nuclei system into a coherent population trapping 
state by optical excitation induced nuclear spin diffusion can be described in terms of Levy flights, 
in close analogy with sub-recoil laser cooling of atoms. The large difference in electronic and nuclear 
time scales simultaneously allow for a measurement of the magnetic field produced by nuclear spins, 
making it possible to turn the lasers that cause the anomalous spin diffusion process off when the 
strength of the resonance fluorescence reveals that the nuclear spins are in the desired narrow state. 



The phenomenon of coherent population trapping 
(CPT) in three-level emitters [T] is at the heart of a num- 
ber of key advances in quantum optics, such as sub-recoil 
cooling of atoms [2] and slow-light propagation [3HS]- In 
these experiments, optical excitation from two low en- 
ergy (spin) states to a common optically excited state 
vanishes due to a quantum interference effect, leading to 
the formation of a dark resonance whenever the two driv- 
ing laser fields satisfy the two-photon resonance condi- 
tion. The fundamental limit on how well quantum inter- 
ference eliminates optical absorption is provided by the 
decoherence rate of the two low-energy spin states. Typ- 
ically, this decoherence rate is assumed to be induced by 
a reservoir which could be treated using the usual Born- 
Markov approximation, implying that the reservoir has 
a short correlation time and its density operator is not 
influenced by the interactions. 

Unlike their atomic counterparts, solid-state spins are 
in general subject to non-Markovian dephasing [5H5] due 
to their coupling to reservoirs with long correlation times. 
In particular, hyperfine coupling to nuclear spins consti- 
tutes the most important source of decoherence for spin 
qubits. ft has been proposed that polarizing or cooling 
nuclear spins could alleviate this decoherence process [6] , 
which prompted theoretical [HI HO] as well as experimen- 
tal efforts aimed at narrowing down the Overhauser field 
distribution [TTlll3| . These schemes could be considered 
as a form of reservoir engineering; remarkably, recent ex- 
periments showed that the substantial manipulation of 
the nuclear spins (reservoir) could be achieved by using 
the electron spin (system) itself [L2rll7| . 

In this Letter, we show that CPT in the spin states 
of a solid-state emitter could be used to prepare a nu- 
clear spin environment in states with a near-deterministic 
Overhauser field. The preparation of an ultra-narrow nu- 
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clear spin distribution is achieved by optical excitation 
induced anomalous diffusion processes [2J. As a conse- 
quence of the anomalous diffusion, the coupled electron- 
nuclei system dynamically switches back and forth be- 
tween a trapped regime where nuclear spin diffusion slows 
down drastically due to the formation of a dark state, and 
a non-trapped regime where optical excitation leads to 
fast diffusion [TS] . We find that for a range of system pa- 
rameters, the coupled system finds the dark state via this 
diffusion process and then remains trapped in it for long 
times, ensuring a narrow nuclear spin distribution with a 
standard deviation that is close to the single-spin limit. 
An additional remarkable feature of the scheme that we 
analyze is the possibility of using resonantly scattered 
photons to measure which regime the coupled system is 
in at a given time [IS]; turning the laser fields off after 
determining the coupled system to be in the dark state 
can then be used to further narrow down the Overhauser 
field distribution to the sub-single-spin regime. Such a 
feedback mechanism is enabled by the large difference 
between the time scales corresponding to electronic light 
scattering and nuclear spin-flip processes. 

Nuclear-spin selective coherent population 
trapping 

We consider a solid-state emitter where the two ground 
electronic spin states, denoted by \f x ) and \l x ), are cou- 
pled by two laser fields to a common optically excited 
state \t) (Fig. la). The laser field with frequency uj p (uj c ) 
that couples the |t x ) — \t) (|J, X ) — \t)) transition with Rabi 
frequency Q p (fi c ) is referred to as the probe (coupling) 
field. The state \t) decays in turn via spontaneous emis- 
sion back to the two ground spin states with an equal 
rate = T t i — Y/2. Denoting the Zeeman energy of 
the electron spin due to the external field B x with lo x 
and the energy of the optically excited state with u t , we 
express the bare optical detunings relevant for the CPT 
system as Aw p = uj t — uj v and Aw c = w ( - w x - w c . 
In the absence of any spin interactions or decoherence, 
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laser fields satisfying the two-photon resonance condition 
(5 = Aujp — Alo c — 0) pump the electron spin into the 

dark state \D) = |t„) - \± x ) , winch 

is decoupled from optical excitation. When Acj c = 
and Q, p ,Q c -C T, the absorption lineshape of the emit- 
ter appears as a Lorentzian with a quantum interference 
induced transparency dip in the center, with a width 

<5zw s ~ (n? + n?)/r « r. 

In practice, the electronic spin states of most solid- 
state emitters are mutually coupled via hyperfine inter- 
action with a nuclear spin ensemble consisting of N nuclei 

N ( 1 \ 

#hy P = 9^2 9i \I> X + -(/>_ + r_a + )\ . (1) 

i ^ ' 

Here, <7i defines the normalized hyperfine coupling con- 
stant between the emitter electron and the i th nucleus 
= !)• I n this convention g = An/J^di quantifies 
the collective hyperfine coupling strength, with Ah de- 
noting the hyperfine interaction constant of the material. 
<7fc and P a (a = +,—,£&) are the electronic and nuclear 
spin operators, respectively; er + = |tx) (ix\- 

Our analysis of CPT in the presence of hyperfine in- 
teractions with a nuclear spin reservoir starts with the 
master equation, obtained by eliminating the radiation 
field reservoir using a Born-Markov approximation: 

r 

p=-{t S ® Pit - {\t)(t\,p} + ) - %[H + #lascr + H hyp ,p], 

(2) 

where ptt — (t\p\t) ac t s on the Hilbert space of nuclear 
spins and Is = Its) (tx| + \Xx) (ix\- We assume that in 
the absence of optical excitation, the electron spin is well 
isolated from all reservoirs other than the nuclear spins 
|20j . and spin-flip co-tunneling or phonon emission rates 
are negligible within the timescales of interest. 

In the limit of a large external field (lu x ^> g), the di- 
rect electron- nuclei flip-flop processes I+er_ + I-<r+ (col- 
lective spin operators are defined as I a = J2i9i^a) are 
strongly suppressed due to the large mismatch in the elec- 
tronic and nuclear Zeeman splitting. In contrast, optical 
excitation does allow for energy conservation in an opti- 
cally assisted electron-nuclear spin-flip process. We take 
the higher order processes into account by applying a 
Schrieffer-Wolff transformation to eliminate the direct 
hyperfine flip-flop interaction. The master equation then 
reads 

r 

P = ^ (IS ® Ptt ~ {\t) (t\,p}+) - i[H + fflascr + ^hyp, p] 

+ e 2 ^l s ® D(p tt ) (3) 
=C (p) + e 2 £ 1 (p u ), 
where the new term containing 

D(p) = I+pI_ + I_pl + - + I-I+,p}+ (4) 



describes an optically induced random nuclear diffusion 
process caused by the optically assisted hyperfine flip- 
flop processes, that are lowest order in the parameter 
e = g/(2uj x ). In Eq. Q wc have neglected terms cx e 2 
that only affect the electron evolution [2"T] . 

After the Schrieffer-Wolff transformation the Hamilto- 
nian relevant for electron spin dynamics (to highest order 
in e) is H spin = Hh. yp + Sa x = ga x (I x + eI + I-+5/g). The 
electron experiences an effective magnetic field, which is 
composed of the two-photon detuning <5, as well as a con- 
tribution originating from the nuclei, which we refer to 
as the generalized Overhauser field I x = I x + el + l-. For 
a given laser detuning <5 each eigenvector of the gener- 
alized Overhauser field I x |A) = A|A) corresponds to a 
steady state p\ — p e {\) <£> |A) (A| of the unperturbed evo- 
lution £o{p\) = 0. Here, p e (X) is given as the solution 
of the optical Bloch equations (OBE) found after projec- 
tion of the unperturbed master equation on the respec- 
tive nuclear state |A); in the OBE i5 e ff = <?A + 6 gives the 
effective two-photon detuning that determines the CPT 
condition. The lifetime of such quasi-steady states p\ 
under the full dynamics of Eq. ([3| is determined by hy- 
perfine assisted scattering events, which are described by 
the term C\. The corresponding nuclear spin flip rate 
in positive (negative) direction D + (D~ ) can be directly 
deduced from Eq. 0: D ± = e 2 §/4(A) (I T I±), where 
Ptt(X) = (t\p e (X) \t) is the population in state |t). Each 
nuclear spin flip event of this kind changes (I x ) by a value 
of order g t . 

For nuclear states with gX = —d the system is in 
two-photon resonance and the electronic system is trans- 
parent such that p\ t = 0: as a consequence, the nu- 
clear spin diffusion vanishes and the system is trapped 
in a dark state. Since the generalized Overhauser field 
in an electronic-nuclear dark state is locked to a fixed 
value, its variance will be strongly reduced (nuclear state 
narrowing) suppressing hyperfine-induced electron spin 
decoherence. Strikingly, by narrowing the generalized 
Overhauser field, even electron-mediated nuclear spin dif- 
fusion is suppressed, thus eliminating the second order 
contribution to hyperfine-induced electron spin decoher- 
ence as well. For all nuclear states satisfying gX s» —6, 
the excited electronic state population will remain small 
(pt t oc S 2 S ), ensuring that the spin diffusion rate will re- 
main vanishingly small: we refer to this subspace as the 
trapping region. 

In contrast, nuclear states with gX ^ —6 render the 
electron optically active and the generalized Overhauser 
field experiences random diffusion (recycling region). To 
illustrate the dynamics allowing the nuclei to move from 
the recycling to the trapping region, we consider an elec- 
tron that is optically excited to state \t): as it decays, it 
could induce a nuclear spin flip event (with probability 
~ e 2 ) in either direction. Through successive spin-flip 
events, the nuclear reservoir probes different spin con- 
figurations with distinct generalized Overhauser shifts. 
When the diffusion allows the nuclei to reach a configu- 
ration that yields <5 e ff ~ 0, the electron becomes trapped 
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in the dark state; further optical excitation is then inhib- 
ited and nuclear spin flips are strongly suppressed. 

Owing to the quasi-continuous nature of the gener- 
alized Overhauser field spectrum, the dark-state condi- 
tion 5 c ff = can be satisfied for a wide range of ini- 
tial detunings Au p . This leads to a drastic change in 
the CPT signature in absorption spectroscopy: instead 
of exhibiting a narrow transparency dip at (bare) two- 
photon- resonance (S — 0), the coupled electron- nuclei 
system displays a broad transparency window. 

The operator valued correction I x to the two-photon 
detuning S and the optically induced diffusive dynamics 
of I x described by the second line of Eq. ^ arc at the 
heart of the nuclear-spin cooling scheme we analyze in 
this work. The predictions we outlined hold in general 
for any nuclear operator I x with a sufficiently large den- 
sity of states around gX + S — 0. In the Methods we show 
that this requirement is fulfilled for the generalized Over- 
hauser field I x = I x + eljf-I- and that its properties are 
very similar to those of I x for the parameters we consider. 
Therefore, for the sake of simplicity, we will proceed by 
neglecting the e correction. As a further simplification we 
will constrain our analysis to nuclear spin 1/2 systems. 
While our results apply to a broad class of solid-state 
emitters, ranging from various types of quantum dots to 
NV centers, we will focus primarily on a single electron 
charged quantum dot (QD) where the optically excited 
state is a trion state consisting of an electron singlet and 
a valence-band hole (Fig. la) [H [22123]. For most QD 
systems, the assumptions we stated earlier are realized 
in Voigt geometry where B x is applied perpendicular to 
the growth direction. 

Semiclassical analysis 

We first consider the semiclassical limit to numeri- 
cally confirm the principal striking features of the cou- 
pled electron-nuclei system - altered CPT signatures 
and the drastic nuclear state narrowing - for inhomo- 
geneous electron-nuclear coupling. To obtain a semiclas- 
sical description of the coupled electron-nuclei dynamics, 
we start by assuming that the electron (p e ) and the nu- 
clear (p n ) spins remain unentangled throughout the sys- 
tem evolution (p = p e <g> p n ). Since, as discussed earlier, 
the electron dynamics takes place on a timescale that is 
faster by a factor e~ 2 ^> 1 than the nuclear dynamics, it 
is justified to solve the OBE in steady state to determine 
the trion population pf t for a given nuclear spin configu- 
ration (and the associated effective magnetic field). 

To describe the nuclear spin dynamics semiclassically, 
we assume that the nuclear density operator p n is diago- 
nal in the basis of individual nuclear spin eigenstates. 
This assumption is justified for QDs in which either 
strongly inhomogeneous hyperfine coupling or inhomo- 
geneous quadrupolar fields lead to large variations in the 
splitting of the nuclear spin states; when this is the case, 
the nuclear superposition states will effectively dephase, 
justifying the assumption of a diagonal density operator. 
In this limit, the master equation Eq. (|3J) reduces to rate 



equations which can be numerically solved using Monte 
Carlo techniques (see Methods). 

Figure 2 shows the result of the Monte Carlo simula- 
tions of the coupled electron-nuclei evolution. To obtain 
the probe field absorption lineshape as well as the Over- 
hauser field variance, we assume that for each probe field 
detuning, we start out from a completely mixed p n , take 
Alu c = and evolve the coupled system to its steady 
state for a range of probe laser detunings. We find that 
the transparency window that has a width of ~ 0.1 2r 
(~ 0.48r) for n c = n p = 0.2T (O c = O p = 0.4T) in the 
absence of hyperfine coupling (Fig. 2a, red dashed curve) 
is drastically broadened and assumes a width 5vt Iana > T 
(Fig. 2a, solid curves). This dragging of the dark reso- 
nance effect is in contrast to Faraday geometry experi- 
ments where nuclear spin polarization ensures that the 
applied laser field remains locked to a detuning that en- 
sures maximal absorption [13] . Concurrently, the Over- 
hauser field distribution is narrowed dramatically from 
its value in the absence of optical excitation (Fig. 2b, 
black dashed line) such that its standard deviation a of 
is smaller than the change induced by flipping one nu- 
clear spin of the most weakly coupled class (Fig. 2b, solid 
curves). These simulations show all the striking features 
that are a consequence of the optically induced nuclear 
spin diffusion [second line of Eq. ([3|] which leads to a 
uni-directional evolution into the electronic-nuclear dark 
state pn = \D)(D\<S> P£), where p 1 ^ is a nuclear spin den- 
sity operator that yields 5 e g = 0. 

We remark that the narrowing of the Overhauser field 
distribution could be measured by using the same two- 
laser set-up and scanning the probe laser on timescales 
short compared to those required to polarize the nuclear 
spins, thanks to the large separation between the elec- 
tronic and nuclear dynamical timescales. Figure 2c shows 
the simulation of the absorption lineshape obtained by 
such probe laser scans. Starting out with 100 random nu- 
clear spin configurations, we first let the coupled electron- 
nuclei system evolve to its steady-state under initial laser 
detunings Aw c = and Aw p = — 0.2r. We then scan 
the probe laser in either direction to obtain the absorp- 
tion lineshape (solid blue curve) that directly reveals in- 
formation about the narrowing of the Overhauser field 
distribution. If we repeat the numerical experiment by 
assuming that the laser fields were initially completely 
off resonance, we find that the absorption lineshape is 
nearly Lorentzian (dashed green curve). 

Quantum model of the nuclear spin dynamics 

Next, we study the homogeneous coupling limit using 
a full quantum treatment. To capture the full quantum 
dynamics, we derive a master equation which depends 
only on nuclear degrees of freedom, allowing for both 
an analytical steady state solution and the comparison 
between the quantum and the semiclassical limit. To this 
end, we assume the homogeneous limit (<?i = 1/y/N). 
First, we eliminate the state \t) in the limit 17 p , Sl c <C 
F, giving a master equation involving the nuclear and 
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electronic spins only. We also assume 
which ensures that the relevant electron spin states in 
the rotating frame are \D) = | l z ) and \B) — \ f z ), and 
choose 6 = 0, for simplicity. In the interaction picture we 
then obtain from Eq. ^ the reduced master equation 

p=r cff (*-P<? z + -lw z + v-,p}+) (5) 

where a z ± are the electron spin matrices in the z-basis 

and r o ff = (r/2) 2 +(gj /2) 2 § * s an °P era tor valued effec- 
tive (electron) spin decay rate. The last line of Eq. ^ 
describes the nuclear spin diffusion determined by the 
nuclear operator p tztz = (f 2 \p\ tz), cf. Eq. Q [25]. 

In order to eliminate the electronic degrees of free- 
dom from Eq. ^ we once again make use of the fact 
that on the timescales of the electron evolution, the nu- 
clear field can be considered as quasi-static and hence 
the electron settles quickly (on nuclear timescales) to its 
interim steady state. We find that on this coarse grained 
timescale = |[1 - (j2^Jj) 2 ]/On, with |A off | 2 = 

r 2 ff + (gl x ) 2 . Using this relation, the electron spin can 
be eliminated from Eq. yielding 

p'» =Tr s (p) = D(T nuc p n ), (6) 

where we defined the nuclear spin flip rate T nuc = e 2 [1 — 
( | X° " | ) 2 ]r c ff ■ Note that since T nuc vanishes for all zero 
eigenstates of I x , Eq. ^ implies, in accordance with the 
considerations above, that every state in the kernel of the 
collective nuclear spin operator I x - i.e. a state of vanish- 
ing Overhauser field - is a steady state of the dynamics. 
We plot the nuclear spin flip rate in Fig. 3a. 

Equation ^ can be used to directly compare the quan- 
tum mechanical and semiclassical diffusion rates in the 
homogeneous limit (see Methods). Surprisingly, the two 
opposite regimes of semiclassical and quantum mechan- 
ical description show both qualitative (evolution can be 
fully characterized by rate equations) and quantitative 
(for the relevant states the calculated rates are compa- 
rable) agreement (cf. Fig. 3a); this result is particularly 
interesting since we would expect the semiclassical de- 
scription to fail in the homogeneous limit. 

In order to calculate the achievable Overhauser field 
standard deviation oof we numerically compute the ex- 
act steady state solution of master equation ^ for ho- 
mogeneously coupled nuclei. To this end, we explicitly 
consider all orders of the hyperfine interaction including 
processes that result in a (small) finite decay rate out of 
the dark state [23] ■ Figure 3b shows a of as a function 
of Q, where we find that <jof decreases with decreasing 
fi until it reaches a minimum of a of — 2 (<t of — 0.7) 
for £1 ~ 0.2r and an electron spin decoherence rate of 
T-2 1 = 100s- 1 (17 1 = 0). This result can be under- 
stood by recalling that the width of the transparency dip 



in CPT scales as fi 2 /r, implying that the range of Over- 
hauser field values yielding transparency can be narrowed 
simply by reducing ft. For fl < 0.2T, we find that a of 
increases rapidly; for such small values of O, the cou- 
pled electron-nuclei system spends substantial amount of 
time outside the narrow transparency region, leading to 
the observed increase in steady state value of a of- As 
we will argue below, this increase does not constitute a 
fundamental limitation for the attainable narrowing and 
values oof < 1 are possible by using a feedback mecha- 
nism. Clearly though, such a remarkable level of narrow- 
ing could only be observed if it is achieved on timescales 
short compared to those imposed by electron spin deco- 
herence and optical excitation independent nuclear spin 
decay processes. We now turn to the question how quickly 
the nuclear spins reach this narrowed state. 

Evolution of the nuclear spins as Levy flights 

There is close analogy between the problem of CPT 
in the presence of inhomogeneous hyperfine interactions 
with a slow nuclear spin ensemble and that of a one- 
dimensional velocity selective CPT [3] [T8]; the role of 
atomic momentum in the latter case is assumed by the 
nuclear Overhauser field I x in the present problem. Just 
like the atomic momentum along the direction of interest 
could change by any value up to the full recoil momentum 
upon light scattering, the Overhauser field could change 
by any value, thanks to an inhomogeneous distribution 
of hyperfine interaction constants < Qi < gl nax . The 
two models differ in two important aspects: first, there 
is a maximum value of the Overhauser field (gl x ) = Ah 
given by full polarization of the nuclei, and second, only 
a small fraction e 2 C 1 of light scattering events give rise 
to a change in the nuclear spin configuration. 

It is known from the velocity selective CPT problem 
that the timescales for sub-recoil cooling of the atomic 
momentum distribution could be understood using Levy 
flight analysis [15]. We apply this method to determine 
the timescale over which we expect the nuclear spins to 
reach a configuration with a mean Overhauser field that 
is smaller than a prescribed value (trapping region). Con- 
sider the random walk of the Overhauser field in time for 
AiOp = Aw c = 0, fl p = Q c = O: this random walk is 
characterized by periods of diffusion followed by long in- 
tervals where the Overhauser field is restricted (trapped) 
at a value close to (I x ) = 0. The duration of the longest 
trapping interval is typically on the order of the interac- 
tion time - a signature of Levy statistics. The probability 
distribution functions P(t) for the trapping time t, and 
P(t) for the recyling time t for which the Overhauser 
field diffuses within the recycling region before returning 
to the trap characterize the Levy flights [18] . 

We are mainly interested in the time required for an 
Overhauser field initially in the recycling region to diffuse 
to a value that is within a prescribed interval that defines 
the trapping region. Once again, this simplification is a 
consequence of the fact that a drop in scattered light in- 
tensity could reveal whether or not the nuclear spin dis- 
tribution has the prescribed value on timescales smaller 
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by e 2 than those needed to flip another nuclear spin. A 
feedback mechanism could therefore ensure that laser ex- 
citation is turned off and the desired/attained a of is pre- 
served. To simplify the estimation of the trapping time, 
we consider a limiting case where f2 2 /V < Ah/N; i.e. a 
typical single nuclear spin flip will take the system out of 
the transparency window. 

Around the transparency point, the rate at which nu- 
clear spins flip is given by D ± oc (I x ) 2 ■ This dependence 
yields P(t) oc t _5 , which in turn leads to infinite aver- 
age trapping times 18J. If we assume that the width 
of the recycling region is determined by ^= ~ then 
the light scattering rate outside the transparency window 
could be taken to be constant with value 2 /r [27]. The 
nuclear spin flip rate in the recycling region is then given 
by r^ 1 « e 2 tf/T. 

In this simplified model the random walk is confined 
and unbiased so in the limit of many nuclear spin flips, 
the number of steps required to return to the trap is 

given by (M) = a ^ff. = y/N . Since the 

time for a single spin flip is taken to be independent of 
the Overhauser field, the time to return to the trap is 
given by 



(t) = (M)r 



a h /vn r 1 
n 2 /r fi 2 ? 



A~h~T 2 



(7) 



For uj x rj Ah this expression simplifies to (t) « /Ah- 
Given the strong N dependence of (t) corresponding 
to the timescale needed to establish a of ~ 1, it is im- 
portant to consider nuclear spin dynamics arising from 
optical-excitation-independent nuclear spin diffusion or 
decay processes, as well as the electron spin decoherence. 
The ultimate limit for the latter is due to spin-orbit me- 
diated spin-flip phonon emission with a rate ~ 10 T for 
uj x ss Ah [28]; as seen in Fig. 3b the resulting increase 
in a of is a factor ~ 3 as compared to the case with no 
electron spin decay. Physical processes leading to nuclear 
spin diffusion in the dark state include (a) nuclear spin 
diffusion mediated by exchange coupling of the QD elec- 
tron spin to a degenerate electron gas or by phonon emis- 
sion/absorption 8J, (b) electric field fluctuations in the 
QD environment leading to spatial shifts in the electron 
wave-function, (c) nuclear quadrupolar fields with axes 
not parallel to B x . If we denote the optical excitation 
independent single nuclear spin diffusion rate that can 
arise from any of these mechanisms with j n and assume 
that N~f n <C ri 2 /r, then we could write the steady-state 
standard deviation of the Overhauser field as 



oof 



A 



H 



(t) + (t) vn (t) + (i) 



(8) 



where the average time spent in the trapping region 
(t) = (Njn)^ 1 and the effective width of the trap 
8 = e~ l Vtyf N'Yn/T. The optimal condition for measure- 
ment independent Overhauser field narrowing is obtained 
when the two contributions to oqf are comparable. As 



we have argued earlier, the use of feedback from the res- 
onance fluorescence intensity should allow for reaching 
oof ~ S. 

Prospects for experimental realization 

We have seen that the strength 7„ of optical excitation 
independent nuclear spin diffusion processes determines 
the degree of attainable Overhauser field narrowing. In 
this context, we remark that experimental observations 
reported by the Bayer group |29] , obtained by driving an 
ensemble of single-electron charged QDs using periodic 
ultra-short optical pulses in the Voigt geometry, demon- 
strated that optically prepared nuclear spin states could 
survive for ~ 10 minutes |30j . Such long nuclear spin 
lifetimes in principle allow for reaching a of ~ 1 using 
the proposed CPT scheme. We also note that the ba- 
sic signatures of CPT have been observed in both single 
electron [35] and hole [53] charged QDs. 

Even though we have concentrated on nuclear spin dif- 
fusion associated with the ground-state hyperfine cou- 
pling, the conclusions of our work remain unchanged if 
the solid-state emitter has hyperfine coupling leading to 
nuclear spin diffusion in the optically excited state. This 
would be the case for example in QDs with vanishing 
heavy-light hole mixing leading to near-resonant hole- 
mcdiatcd nuclear spin-flips in the excited state due to 
the dominant S^I Z term in the hole-hyperfine interac- 
tion Hamiltonian. 

While prior experimental results on pulsed excitation 
of an ensemble of QDs strongly suggest the feasibility of 
our proposal in self-assembled QDs, we expect our find- 
ings to be relevant for a wider range of solid-state emit- 
ters. Of particular interest is nitrogen- vacancy (NV) cen- 
ters in diamond where CPT has also been previously ob- 
served |31] . The small number of nuclear spins coupled to 
the optically excited spin in the case of NV-centers should 
make it possible to reduce the time needed for the system 
to find the dark state drastically. A principal difference 
with respect to the large N limit we analyzed is the fact 
that only a small set of optical detunings will allow the 
NV system to find a dark state. Finally, extensions to 
other solid state systems such as superconducting qubits 
may be possible [32] . 
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METHODS 

Hamiltonian of the laser driven coupled 
electron-nuclei system. 

The Hamiltonian of the solid-state CPT system we 
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consider is He pt = Hq + -ffiascr + -^hyp with 

H = oj x a u + uj t a tt , (9) 



H lasei = [lVn e ~ l " pt + "c^e- iWct + h.c 



(10) 



Here, f (I) denotes fx (Ix)- The definition of Hh yp as well 
as the quantities appearing in Hq and -ffiascr are given in 
the main text. 

Rate equation description of nuclear spin dy- 
namics. 

The semiclassical limit can be derived from the mas- 
ter equation by replacing the collective spin decay by 
independent decay of individual spins. This is read- 
ily accomplished by making the substitution I + p n I- = 
Y,t ] 9i9jI l +P n I+ Y,i9il+P n l\ [and correspondingly 
for other terms in Eq.Q]. We coarse grain the nu- 
clear motion with regard to the electron dynamics and 
from the new master equation we obtain a rate equation. 
We introduce a shell model of the QD with M differ- 
ent classes of nuclear spins (Fig. lb) ; the nuclei in class 
{y) have identical g v and their net spin polarization is 
m v = \{N+-N~) = (E«,4), where N+ (N~) denote 
the total number of up (down) spins in class (v). The 
derived rate equation for the joint probabilities P({m,,}) 
associated with the nuclear spin configuration {m^,} is 
given by 



at 



M 



= E 7 '^})^ ({%})!+({%» 

V 

M 

v 
M 



where Tj_({m fl }) 



"i ' vf^) 2 !/ 9 *^!^}) are the rates at 



which nuclear spins of the i^th class are flipped if the nu- 
clear spin polarizations in each class are given by {m^}. 
{%} ({"v}) denotes the nuclear spin configuration that 
differs from the configuration {m^} only in the vth class, 
with polarization m v — 1 (m„ + 1). 

We numerically simulate the evolution of the nuclear 
spins with a Monte Carlo method. We assume in our 
numerical simulations that the QD contains 100 nuclear 
spins. We group these spins into five concentric shells 
(AI = 5) with different hyperfine coupling constants that 
are determined by the 3D Gaussian electronic envelope 
function (Fig. lb). The coupling constants ggt for these 
shells are 0.0934r, 0.0828r, 0.0678r, 0.0513r, 0.0358r 
and the corresponding total numbers of nuclear spins in 
each shell are 2, 8, 16, 28, 46. The coupling constants 



are chosen to ensure that the standard deviation of the 
Overhauser field seen by the QD electron for nuclei in a 
completely mixed state satisfies <Jof(p) = j- We do not 
keep track of the exact configuration within each class 
(i/) of nuclear spins and assume that any configuration 
of spins leading to the same m M is equally likely and that 
the nuclear spin distribution in each shell is independent 
of the other shells. 

Nuclear spin flip rates for homogeneous hyper- 
fine coupling 

For any eigenstate \m) of I x with I x \m) = m/^/N \m) 
(r nU c | fn) = F™ c |m)) the nuclear spin flip rate in nega- 
tive (positive) direction is given by D~ = (/ + /_r nuc ) m 
(D + = (7_/ + r nuc ) r „). In the semiclassical limit under 
the assumption (P + I J _) = (i ^ j) the rates are sim- 
ply given as Df c = (^ ± ^)r™ c . For the quantum 
description the characterization via the spin projection 
quantum number m is not sufficient; the rates also de- 
pend on the symmetry of the nuclear state, quantified 
by the total spin J e {0, N/2}. For a Dicke state 
\J,m) the rates are given as DJ m = (I±I T T nuc )j^ m = 

^[J(J + 1) — m(m =F l)]r™ c . For the statistically rele- 
vant J-subspaces [J = 0(VN)] we find that D^ m are in 
good agreement with those obtained in the semiclassical 
limit for small m values. 

Generalized Overhauser field 

For the main part of the nuclear Hilbert space - namely 
the domain where the operator I x + S/g is large (recy- 
cling region) - the e correction to the Overhauser field 
represents a negligible perturbation to the hyperfine in- 
teraction. However, in the domain of small eigenvalues of 
I x + S/g (trapping region) this perturbative picture is not 
trivially justified. However, for the diffusive dynamics we 
are mainly interested in the number of eigenstates in a 
1 region L = (—5 — r),—6 + ri)/g around —S/g. For I x the 
number of eigenstates with small eigenvalues is very large 
(exponential in the number of spins), which is favorable 
for our scheme. For the generalized Overhauser field it is 
given by D{ V , e) = JZ^ dE Tr(<5(£ - I x - eJ+J_)) and 
deviates from the number of unperturbed eigenstates: 



(11) 



71 = 1 



The sum ^ri=i can be upper bounded by ~ D(r], 0) 
for large trapping regions r/ 3> e, i.e. the number of eigen- 
states changes at most by a factor of order 1. Numerical 
calculations for N = 10 4 inhomogeneously coupled spins 
show that even for r\ ~ e the number of states in both 
the perturbed and unperturbed case differ only by a few 
percent. 
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Figure Captions 

Figure 1 The energy level diagram of a solid-state 
emitter, (a) The electron spin state | fx) (I ts)) is 
resonantly coupled to a trion state with an x (y) polar- 
ized laser field with Rabi frequency ft p (ft c ). In Voigt 
geometry the oscillator strengths of the two transitions 
are identical, leading to spontaneous emission rates with 
equal strength T/2. Optical excitation allows for energy 
conserving hyperfine flip-flop transitions that result in 
nuclear spin diffusion; these second-order processes are 
depicted using dashed curves, (b) The confined electron 
wave-function leads to inhomogeneous hyperfine coupling 
with the nuclei. In the simulations, we assume that the 
dot can be described as consisting of 5 different classes of 
nuclei. All nuclei within a class have identical hyperfine 
coupling, with strength determined by the electron wave 
function. 

Figure 2 Nuclear spin state selective coherent 
population trapping, (a) The absorption lineshape in 
the presence of hyperfine interactions with quantum dot 
nuclei for Rabi frequencies ft p — ft c — 0.2 V (red) and 
ftp = ft c = 0.4 r (blue): in stark contrast to the stan- 
dard coherent population trapping profile (dashed lines), 
the dark resonance is drastically broadened (solid lines). 
The broadening of the dark resonance is a consequence of 
the fact that optical excitation induced nuclear spin dif- 
fusion allows the coupled electron-nuclei system to find 
a Overhauser field configuration that satisfies the dark 
state condition for a broad range of initial laser detun- 
ings. (b) The standard deviation of the Overhauser field 
a of for ft p = ft c = 0.2 T (red line) and ft p — ft c = 0.4 T 
(blue line)) is reduced to the level below that of a single 
nuclear spin flip (green dashed line). The dashed black 
line shows the standard deviation in the absence of laser 
drive, (c) The final Overhauser field distribution could be 
determined by using a fast scan of the probe laser across 
the optical resonance. The solid (blue) curve shows the 
absorption lineshape obtained when the system is ini- 
tially prepared in a dark state by setting Aoj p = — 0.2F: 
the width and the depth of the dark resonance reveals 
information about the Overhauser field distribution. If 
the experiment is carried out by starting out in a ran- 
dom nuclear spin state, the observed lineshape is close to 
a Lorentzian (dashed green curve). 

Figure 3 Homogeneously coupled electron-nuclei 
system, (a) Nuclear spin diffusion rates depending on 
the nuclear spin projection m assuming homogeneous 
coupling (cf. Eq. (J6|). Parameters are N — 4 x I0 4 , 
r = I GHz, A = ui x = 100 ^cV and ft = 0.1 GHz. 
The quantum mechanical rates are taken for the subspace 
,/ = y/N/2. (b) The dependence of the steady state value 
of the nuclear Overhauser field standard deviation as a 
function of ft, calculated using the fully quantum model. 
The increase in standard deviation for ft„ = ft c < 0.2r is 



a consequence of the fact that the coupled electron-nuclei 
system spends a substantial amount of time outside the 
transparency region. The inset shows the steady state 
population of I x eigenstates for ft p = ft c = 0.02I\ De- 
spite the finite standard deviation a ps 9 the system is 
strongly peaked around m = 0. 
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